Superconducting phase diagram and FFLO signature in A-(BETS)2GaCl4 from rf 

penetration depth measurements 



William A. Coniglio, Laurel E. Winter, Kyuil Cho, and C. C. Agosta 
Department of Physics, Clark University, Worcester, Massachusetts 01610 

B. Fravcl and L. K. Montgomery 
Department of Chemistry, Indiana University, Bloomington, IN 4^405 
(Dated: May 13, 2011) 

We report the phase diagram of A-(BETS)2GaCl4 from rf penetration depth measmements with a 
tunnel diode oscillator in a pulsed magnetic field. We examined four samples with 1100 field sweeps 
in a range of angles with the magnetic field parallel and perpendicular to the conducting planes. In 
the parallel direction, Hc2 appears to include a tricritical point at 1.6 K and 10 T with a phase line 
that increases to 11 T as the temperature is decreased to 500 mK. The second phase line forms 
a clearly defined high field low temperature region satisfying several of the conditions of the Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) state. We show remarkably good fits of IIc2 to WHH in the 
reentrant a > 1, Xso = regime. We also note a sharp angle dependence of the phase diagram about 
the field parallel orientation that characterizes Pauli paramagnetic limiting and further supports the 
possibility of FFLO behavior. Unrelated to the FFLO study, at fields and temperatures below IIc2 
and Tc, we find rich structure in the penetration depth data that we attribute to impurities at the 
surface altering the superconducting properties while maintaining the same crystallographic axes as 

Hc2- 



I. INTRODUCTION 



Organic conductors provide an excellent platform for researching unconventional superconductivity because they 
are clean low-dimensional materials whose critical temperatures and fields are accessible, yet their superconducting 
mechanisms arc not yet fully understood. By developing the Tunnel Diode Oscillator (TDO) penetration depth 
measurement and pulsed field techniques, we aim to discover finer details in the phase diagram of the organics. 

The cation BEDT-TSF (bis(ethylene-dithio)tetraselenafulvalene), usually called BETS, is closely related to the 
commonly studied BEDT-TTF (bis(ethylene-dithio)tetrathiafulvalene), which is shortened to ET, with the innermost 
four sulfur atoms of ET replaced by selenium to form BETS. The A-(BETS)2GaCl4 crystals are needle-shaped and 
crystahze alongside platelets of K-(BETS)2GaCl4, which do not superconduct. Layer spacing for A-(BETS)2GaCl4 is 
18.58 A. [l| A very interesting relationship exists between the title compound and A-(BETS)2FeCl4, which becomes an 
insulator at the title compound's critical ternpcrature. A-(BETS)2FcCl4 also exhibits field-induced superconductivity, 
and may have an FFLO state of its own.0, [11 

In a quasi-2d layered superconductor, the tunneling current between layers is small, but finite, so the layers act 
nearly independently, but maintain long-range coherence. [1] If the vortices are able to fit between the layers, the 
orbital limiting field increases dramatically, and superconductivity persists to much higher fields. 0, Q Barring other 
limits, when the Zeeman energy of a Cooper pair exceeds its binding energy at the Pauli paramagnetic limit, the pair 
will align parallel and the sample will no longer superconduct. 0-13 We have employed the theory of WHIlflol. [Tl| to 
model Hc2{T), taking into account both orbital and Pauli paramagnetic limiting. 

Paramagnetically limited materials present interesting possibilities for novel superconducting behavior. In a sample 
where the orbital limit is significantly larger than the Pauli paramagnetic limit, a large magnetic field may cause the 
ground state wave function to become spatially modulated. In this case, the Cooper pairs acquire a non-zero linear 
momentum. This inhomogeneous superconducting state is known as the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO 
or LOFF) state after the two groups who postulated it independently in 1964. [l2|, [l^ There have been a number of 
reports of FFLO behavior in various samples, although none have been definitive. Even so, the organic superconductors 
remain some of the ripest ground for the search. In the title compound, thermal conductivity measurements show 
a cusp at the appropriate location in the phase diagram. There arc also reports in other organics [Tsl - fl^ . inclu ding 
some using the same technique wc use here. Notably, a-(ET)2NIl4(SCN)4 does not exhibit any FFLO phase. [2C| 
In the heavy fcrmion superconductor CcCoIiib, specific heat measurements [2ll, [2^ suggested FFLO behavior, and 
penetration depth studies [23j followed. Ref. [24 | argues for a competition between ordered states that results in a spin 
density wave. Ref. ^251] suggests that the SDW may coexist with an FFLO state. The actual behavior is still unclear. 
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FIG. 1: Tunnel diode oscillator circuit used in 390 MHz examination of the organic superconductor A-(BETS)2GaCl4. R2 
166 f2, Cbypass = 471 pF, Tip = 100 fi, C = 5 pF. The length of the tuned coax cable is 511 mm. 



II. EXPERIMENTAL DETAILS 

Single-crystal samples of A-(BETS)2GaCl4 were grown using an electrochemical techniqueP, [2^ in 1995-6 and 
stored for 13 years. We will follow the crystallographic convention that the long axis of the needle grows in the 
c-direction, and the best conduction occurs in the a-c plane with layers stacked along the b-axis. Samples A and C 
measured 0.1 mm each along the a- and b-axes and 0.3 mm along the c-axis. Sample D had the same proportions 
in a slightly larger size. Sample B was slightly wider than A and C along the a-axis, electrically different from the 
other three, and does not agree with the literature on A-(BETS)2GaCl4. It appears to superconduct, but at Tc = 4.25 
K and Hc2 = 9.3 T. We measured Sample D (Tc = 5.1 K) in a slightly modified experiment and obtained a limited 
amount of data, all of which was consistent with Samples A and C. 

Penetration depth studies of A-(BETS)2GaCl4 have been difficult in the past due to the small size of the samples. 
We redesigned portions of the tunnel diode oscillator circuit to make effective use of standing waves that appear 
when using high frequencies. We mounted the samples on a 1-axis rotator with 0.1° resolution and ensuring precise 
alignment by pulsing the field up and down at a range of angles 0.23° apart about the parallel orientation. We oriented 
the sample on the rotator to apply the magnetic field along the b- or c-axes (Samples A and C) and along the b- or 
a-axes (Samples B and C). We ran Sample C twice, with fields parallel to the a- and c-axes, and noticed no significant 
differences. Orientation in the a-c plane was set upon sample mounting within a tolerance of ±15°. The magnetic 
field of the oscillator was always along the b-axis, so that currents are excited within the planes of the conductor. 

We performed the experiments in a '^Hc cryostat inside the 40 T medium pulse magnet at Clark University with a 
rise-time of 27 ms. By using either "^He or '^He in the cryostat, we were able to take data from 367 mK to 4.3 K in a 
liquid environment. 

To measure the in-plane penetration depth, we designed a Tunnel Diode Oscillator [27l - [29| shown in Fig. [T] and 
described in (30j . to operate at 390 MHz with the sample contained in a 4-turn 0.5 mm diameter coil. For maximum 
sensitivity, the oscillator circuitry was placed next to the rotator and connected to the sample coil with 18 mm of 
twisted pair. We reduced the amount of heat generated near the sample by moving our biasing and matching resistor 
network above the '^He liquid level and used a length-tuned coaxial cable to connect the oscillator to the matching 
network. We achieved minimal noise in pulsed field despite locating the diode in the region of greatest field by using 
rf techniques such as impedance matching the signal into the 50 O coax leading up the probe, regulating the supply, 
and the use of ferrite chokes around coax terminations throughout. With the LC oscillator free running at 390 MHz, 
our pre-pulse data gave a frequency measurement with a scatter of about ±200 Hz. corresponding to an independent 
0.5 ppm measurement every 13 ps. The noise was about 10^ times smaller than the total frequency shift from the 
superconducting to normal state. 

Our RF receiver used a single-conversion to 5 MHz. The radio environment between 390 and 400 MHz was quiet, 
and conversion images were minimal. The conversion introduces a minus sign into the receiver frequency, and a graph 
of received frequency vs. magnetic field may be interpreted as a penetration depth vs. magnetic field, although it is 
difficult to obtain absolute units for the penetration depth. [3ll,[32| We digitized the signal at 20 MS/s for computational 
frequency demodulation and analysis. [30| 
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FIG. 2; (Color online) Example frequency vs. field at 512 mK with two derivatives shown. At low field, the sample is in the 
superconducting state. As the field increases, Cooper pairs are destroyed, and the London penetration depth Xl increases. 
At 9.5 T (lower arrow), the sample enters the inhomogeneous FFLO state. As the magnetic field increases further, the 
inhomogeneous state can no longer support the magnetic field, and the sample transitions to the normal metallic state, (upper 
arrow) When the penetration depth is larger than S, the metallic skin depth, the flat normal state dominates the RF response. 
Also shown is a fleld sweep at the same temperature from Sample C, which exhibits transitions at identical fields. 



A. Definition of Hc2 



There is inherent ambiguity in the definition of Hc2 because the transition is often broad. We favor identifying the 
critical field using minima in (Pf/dB^, an approach that is well-defined, highly localized, requires little to no human 
judgement, and produces self-consistent results over many trials. We use the exact same criteria (minima in d^f/dB^) 
to identify the Vortex SC-FFLO, FFLO-Normal, and Vortex SC-Normal phase transitions. Inspection of the top plot 
in Fig. [2] clearly show two complete transitions at low temperature in both samples - the lower one from a vortex 
superconductor to FFLO, and the upper one from FFLO to a normal metal. 

In addition, we have found that the critical temperature we determine by an analogous method (minima in cP f /dT"^) 
agrees with the critical fields extrapolated to Tc- Our Tc agrees with a specific heat study [33]. 

When in the superconducting state, the sample excludes RF energy from its center. The excluded volume lowers 
the inductance of the measurement coil, raising the oscillation frequency. [32| In the superconducting state, we measure 
the London penetration depth (Al). In the normal state, the rf exclusion is governed by the metallic skin depth of 
the sample ((5), which depends on frequency and conductivity. From normal state resistivity measurements [ij], we 
expect skin depths from 25 to 80 pm. We measure _ffc2 sX the point where the rapidly increasing penetration depth 
reaches the same size as the skin depth. We can represent that point very precisely by taking the second derivative 
of oscillation frequency as a function of field or temperature. We define Hc2 at the point where the concavity is a 
minimum, signifying a corner between metallic and superconducting behavior. 



III. RESULTS AND DISCUSSION 



We examined each sample in zero magnetic field and with the conducting layers both parallel and perpendicular to 
a magnetic field. Average Tc was 4.93 K, in agreement with (ssl. |33|. From our measurements, we have two reports 
to make. First, we will establish the phase diagram with the field perpendicular to the conducting layers using a 
full superconducting phase diagram measurement on a single sample and argue against the conclusion in [3^ that 
A-(BETS)2GaCl4 exhibits a dimensional crossover at 1.9 K. Second, with the assistance of the WHH model, we will 
use our new penetration depth data with field parallel and nearly parallel to the conducting layers to identify the high 
field low temperature phase as a function of temperature, field, and angle, and affirm the suggestion of |l4l | that the 
phase may be due to the Fulde-Ferrell-Larkin-Ovchinnikov state. We observe that as the angle is rotated away from 
parallel to the planes, the transition between the vortex state and FFLO increases in field, while the introduction of 
vortices destabilizes the FFLO phase and causes its transition to the normal state to reduce in field. We attribute 
the vortex-FFLO behavior to an increase in spin-orbit scattering caused by vortices entering the conducting layers. 
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FIG. 3: H^2{T). Each point represents an independent field sweep, with 53 sweeps total. Points shown are from the upsweep 
only, as vibration introduced extra noise into the downsweep. There was no hysteresis between sweeps. The slight difference 
in critical fields above and below 1.3 K reflects the use of liquid "^He or ^He. Because of the upward concavity above 2.5 K, we 
were unable to realize an acceptable WHH fit in the perpendicular orientation. 



A. Perpendicular phase diagram 

Fig. 13] shows as a function of temperature for a single sample. Comparing our phase diagram to Fig. 9 of [3^ . 
we do not see the abrupt change in slope of Hc2 that they report at 1.9 K. Instead, we see a gradual increase in 
slope at about 3 K. Ref. [s^ suggested that at 1.9 K, the sample undergoes a dimensional crossover from a three- 
dimensional behavior to a quasi-two-dimensional one when the across-plane coherence length becomes smaller than 
the interlayer spacing. While we cannot rule out that explanation, we disagree with their conclusions regarding the 
dimensional crossover. The use of to measure is not easily justified. Whether the sample is 2d or 3d, the 
in-plane size of a vortex depends on ^|| , and their cross-plane extent is the entire sample. Fields aligned parallel to 
the layers should be used to study the anisotropy and dimensionality, as we have done in Fig. [5] and discuss in the 
section on angle dependence of i?c2- 

We find two-dimensional behavior at both 1.35 K and 2.85 K. Since we do sec a gradual change in slope near 3 K, 
we wish to keep open the question about the origin of this behavior. Due to the change in slope, WHH does not fit 
this data. WHH may not generally apply to the organics in the perpendicular orientation because samples are in the 
clean limit, where the mean free path is much greater than coherence length. Using our perpendicular orientation 
phase diagram and extrapolating to zero temperature, we find ^^^g = 106A. 



B. Parallel phase diagram and WHH fit 

The phase at high field and low temperature has been ascribed to the FFLO state. [l3l Our study has produced Hc2 
and Hfflo values with significantly more detail than the previous work. Points from both measurement techniques 



of [iJI and the present work appear in Fig. 01 

We will now compare our data to the theory of Werthamer, Helfand, and Hohenberg[10|, |11|. Starting from fixed 
Tc = 4.9 K and no spin-orbit scattering {Xso = 0), we use the low temperature data to tune Hp = 9.5 T and high 
temperature data to tune 

= o.692r, y-^)^ " "TI " ^^'^ 

From ^-nd Hp, we find a = 3.5 and produce the fit shown in Fig.|4l Remarkably, the fit reproduces the maximum 



Hc2 at 1.6 K, in excellent agreement with both the present work and the interpolated value in |14| . For comparison with 
[35| . using the Dingle temperature from [sj] (whose samples came from the same batch as ours) and our sJ^^q = 106 A, 
we find 3.5 < r < 12. Using only data from [SJl, r might be 10% lower. 

Despite the good fit, we pay careful attention to the intended purpose of the WHH model: First, WHH was 
designed for bulk superconductivity, and no consideration is made for the effects of layers, particularly the orbital 
lock-in effect. We considered the theory of Klemm, Luther, and BeasleyjH, which does examine the effect of layers. 
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FIG. 4: (Color online.) H^2i'^) from 116 magnet sweeps on Samples A and C. The two samples produced data that were 
indistinguishable from one another, and we plot them together here. We see no hysteresis in critical fields between the up and 
down sweeps. Below 1.7 K, a second minimum in the d^f/dH^ representation appears between 9 and 10 Tesla, which forms 
the lower branch on the phase diagram from the vortex superconducting state to the FFLO. Below 1.3 K in one of the samples, 
the Vortex SC-FFLO transition was partly masked by other features in the penetration depth, which we discuss later. Also 
shown are our WHH fit and the transport data from [14| . 



At low temperatures in the parallel orientation, the behavior of KLB and WHH is identical, and WHH require fewer 
parameters. While we eontinue to examine the application of KLB in hopes of modeling the microscopic behavior and 
tilted fields more carefully, the present parallel orientation data appears to be modeled well as a bulk-2d anisotropic 
superconductor with Tc = 4.9 K, a = 3.5, Hp = 9.5 T. Xso ~ 0, and the addition of FFLO behavior at temperatures 
below the reentrant point. In particular, WHH reproduce the maximum in Hc2, which appears to be the tricritical 
point for the FFLO phase. According to MakiQ, that local maximum is also the temperature below which the 
transition becomes first order. Determining the order of a transition from penetration depth data is difficult, and we 
cannot make any strong statements. We instead point out that the present data has a qualitatively different shape 
from CeCoIns, which has a second order transition at the lower field (change in slope of Al in [2^) and a first order 
transition at the higher field (sharp increase in in [1^). 

Second, the model is stated to apply in the dirty limit, which sounds as though it should not apply to the clean 
organic superconductors. With the field aligned exactly in the parallel orientation, cyclotron orbits prevent electrons 
from traveling great distances in most directions without encountering a layer boundary. The presence of layers 
scattering normal electrons causes their mean free path to become short, and the dirty limit becomes an acceptable 
approximation. 

Third, WHH require that spin-orbit scattering is a small portion of the total scattering. Since we must set the 
spin-orbit scattering parameter (A^o) to zero, we have satisfied this limitation of the model. We suspect the orbital 
lock-in effect prevents scattering by removing the vortices from the strongly conducting layers and away from (in direct 
space) the majority of the Cooper pairs. Support for the suppression of spin-orbit scattering in the parallel orientation 
comes from Figs. [5] and [6l where we observe the transition from a vortex superconductor to FFLO (9.7 T in Fig. [5]) 
in the parallel orientation rising as the angle is rotated very slightly away from parallel, and vortices begin to cut 
through the layers. When vortices enter the conducting layers, spin-orbit scattering begins to occur, and the critical 
field bounding the vortex superconducting state rises. The introduction of spin-orbit scattering in the WHH model 
indeed causes an increase in Hc2, even after the reduction in a due to anisotropy is considered through Tinkham's 2d 
equation. Other effects due to the orbital lock-in effect are also possible and are under consideration. WHH should 
not apply when the fields are aligned far from from the parallel orientation, as the mean free path extends into the 
clean limit. The model docs not fit our perpendicular data at all. 

We have attempted preliminary fits to other superconductors with reported FFLO states. While the transition from 
vortex superconductivity to FFLO in the organics seems to fit WHH, CeCoInsjl^ cannot be fitted to a reentrant 
WHH model because the tricritcal point is too low in temperature. The difference suggests that the mechanism for 
the high field, low temperature state in CeCoIns is different from the organics. The competition between ordered 
states explored in CeCoIns [2^ therefore does not discount the possibility of FFLO in the organics. 
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FIG. 5: (Color online.) Hc2 vs. angle from parallel orientation. G-L 3D and Tinkham 2d[33| curves from least-squares fits are 
shown with fit data extending to 90°. The 2d ones are sharply peaked, while the 3D are flat at 0°. The perpendicular end 
of the fits were fixed at measured values of H^. The 2d fits are both visually and statistically better at both temperatures, 
eliminating the possibility of a 2d/3d crossover between them. Parameters at 1.35 K are: 10.44 and 1.96 T, and at 2.85 K, 8.40 
and 0.86 T. We avoid fitting the parallel data, because the fits do not account for the orbital lock-in efi'ect, which appears as a 
nearly discontinuous increase in Hc2 within 2° of parallel. Inset: Full range of data with 2d fit. 
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FIG. 6: (Color online.) Here we illustrate the boundary of FFLO behavior as the temperature is lowered and field angle rotated 
away from parallel to the layers. Points connected by lines all share the same temperature. The dimension in and out of the 
page is the angle of the magnetic field from the layers. The reader may also consider it in 2d, where each group is at the 
X-axis temperature, and points connected by lines show the angular dependence of the critical fields with the apex at exactly 
parallel. The region between the blue and red markers is the FFLO state. At each temperature, the highest field FFLO-Normal 
transition and lowest field Vortex SC-FFLO transition both occur at the same angle (same field sweep, defined as 0°). As the 
angle is rotated in either direction away from 0°, the FFLO-Normal (red triangle) transition lowers in field, while the Vortex 
SC-FFLO one (blue squares) rises to meet it. See text for explanation of the rising behavior due to spin-orbit scattering. The 
angle sweep from Fig. [5] appears here at 1.35 K. 



C. Hc2 vs. angle 

With Sample A, we measured Hc2 vs. angle for a 91° range from slightly before the parallel orientation to the 
perpendicular for temperatures 1.35 and 2.85 K. Resolution of the rotator was 0.12°, and our finest step size was 
0.23°. The portion of both sweeps below 14° is shown in Fig. [5l At both temperatures, even when the immediately 
parallel data is ignored, fitting to the Tinkham 2d form is superior to the G-L 3D one, eliminating the possibility of 
the 2d/3d crossover described in [33|. For small angles about the parallel, we see evidence of the orbital lock-in effect, 
as critical fields jump above the Tinkham prediction below a critical angle. 

In order to ensure the exact parallel alignment of our sample to the magnetic field, we performed field sweeps at 
a range of angles near 0°. Critical fields as a function of angle and temperature appear in Fig. [5] as a 3d plot. The 
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FIG. 7: (Color online.) Our concavity map from Sample A indicates how easily we can identify portions of the phase diagram. 
The color key at right corresponds to the second derivative of oscillation frequency, (fif/dB^. Dark areas trace out the 
superconducting transition, in analogy with Fig. ID White areas are zero concavity. Hc2 is clearly visible as the darkest band 
from bottom-right to top-left. The transition between the vortex state and FFLO branches down from Hc2 at about 1.6 K and 
10 T. The transition-like stripe at lower fields and temperatures is unexpected, and may reflect effects at the surface of the 
sample. 

FFLO state is sideways pyramid-shaped pocket with its point at the tricritical point (1.7 K, 0°, 10.0 T). Of particular 
experimental interest, the critical field between the superconducting vortex state and FFLO rises as the angle is 
rotated away from parallel, even as Hc2 at 11 T falls. 

D. Features below Hc2 

We have noticed features in our data below Hc2 and Tc- We attribute them to impurities at the surface of the 
sample, but their similarity to a superconducting transition interests us. One of them is visible in Fig. [7] from 3.75 K 
to 8.5 T (green online), and also appears in Sample C. A second only occurs in Sample C and matches exactly with 
the only transition visible in Sample B. The newer (and probably cleaner) Sample D has no extra features. We have 
noted similar behavior in three samples of LiFeAs,[33| suggesting an effect of surface reconstruction or oxidation. The 
effect is visible as a deviation of penetration depth from the bottom of the diagonal lines in Fig. 1 of Ref. [s^l • 

Lower features in both the title compound and LiFeAs appear very much the way a superconducting transition 
would, but do not appear to affect the bulk superconductivity in the material. They exhibit no hysteresis between the 
up and down sweep, often have Tc-like signatures in a zero-field temperature sweep, and have an angle dependence 
that follows the 3d G-L or 2d Tinkham forms, depending on the dimensionality of the sample. The crystallographic 
axis is shared with the main transition at higher field, eliminating the possibility of a twinned crystal or a different 
morphology from consideration. 

The coexistence of two forms of superconductivity that appear not to interact but share the same crystal structure 
is puzzling. If the surface of the samples exhibit different superconducting properties, the penetration depth measure- 
ment technique is sensitive to the behavior of one order parameter even if the other is fully superconducting, as the 
susceptibility of the sample will still change. Electrical conductivity, in particular, will not be sensitive to either the 
low field transition or the one between the ordinary superconducting state and FFLO one because the resistance will 
be zero below Hc2- 

IV. CONCLUSION 

We show a phase diagram for A-(BETS)2GaCl4 with 116 field sweeps showing a high field, low temperature phase 
we attribute to the FFLO state. The angular dependence of the phase diagram is consistent with expected FFLO 
behavior. Our great point density allows us to reliably fit the ground state superconducting phase to the model by 
WHH. To our knowledge, this is the first example of experimental data that fits the reentrant WHH phase diagram. 
As no conclusive experiments have been done so far on the FFLO phase in the organics, it is difficult to prove that 
explanation. However, the body of evidence for the phase is growing, and it has been reported from a number of 
different measurement techniques. In addition to the ongoing phase diagram work by several groups, parallel-aligned 
high field studies of specific heat, neutron scattering, or STM, although difficult, would provide valuable confirmation 
or redirection of the research effort so far. 
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